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p ■ ABSTRACT: We present the results of a numerical simulation aimed at understanding 

(-H I the nature of the 'c = 1 barrier' in two dimensional quantum gravity. We study multiple 

Ising models living on dynamical (p graphs and analyse the behaviour of moments of the 

^ ■ graph loop distribution. We notice a universality at work as the average properties of 

c^ I typical graphs from the ensemble are determined only by the central charge. We further 

argue that the qualitative nature of these results can be understood from considering the 

effect of fluctuations about a mean field solution in the Ising sector. 

In recent years considerable progress has been made in understanding the nature of 
some simple systems encorporating two dimensional quantum gravity. Two approaches 
have been followed which are at first sight quite different. The first of these employs 
techniques borrowed from conformal field theory to calculate the spectrum of anomalous 
dimensions in several simple models [1-2]. In the second method a regularisation of the 
continuum functional integrals is made in terms of random triangulations [3-5]. Where 
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it has been possible to solve the models analytically the two approaches have been in 
complete agreement. However, in both cases all attempts to solve the models when the 
matter central charge c becomes greater than unity have failed. The lattice regular isation 
however has the advantage that it allows these 'strongly coupled' regions to be explored 
using computer simulation techniques. 

In this letter we present the results of one such study aimed at examining whether the 
phase transition predicted for c > 1 is a genuine physical phenomena or merely an artefact 
of our current methods of solution [6]. A single Ising model, at criticality, can be shown to 
possess c = 1/2, so a set of P noninteracting systems would (at least naively) correspond 
to a central charge c = P/2. By tuning P from one through three we can search for any 
sign of discontinuous behaviour associated with such a transition. The partition function 
we study can be written as 

Z(/3,P) = 5^e^^(/^'G) 
Gn 



y^ g/3GyCT,crj ^-^^ 



The sum over (j) graphs G^y is restricted to those with A^ nodes and fixed spherical 
topology. The connectivity matrix Gij can be taken to be unity when the sites i and j are 
neighbour on the graph and zero otherwise. 

We sample the space of all graphs G^y by random moves which are dual to direct 
lattice link flips as described in [7] , whilst cluster algorithms are used to provide an efficient 
updating procedure for the Ising sector [8] . On a fixed graph this partition function would 
simply factorise into P copies of a single Ising partition function, but one expects that the 
use of a dynamical lattice induces an effective coupling between different Ising 'species'. 
To this end we collected Monte Carlo data for the standard Ising observables such as the 
specific heat C, susceptibility x, and generalised spin-spin correlation matrix C^^ (r) 
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C"^r) = (a-(0)a/5(r))^ (2) 

The indices a and (3 refer to the species and r measures the geodesic distance on the graph 
between the spins. The latter quantity, in particular, could be expected to show signs of 
this induced coupling by the occurrence of non-vanishing off diagonal entries. However, 
in spite of relatively good statistics (5 x 10 sweeps for a range of (3 from 0.3 to 1.2) we 
were unable (within the statistical errors) to extract any reliable signal that would indicate 
behaviour inconsistent with that of a single Ising model. As a sample of our data we have 
listed in table 1. the effective masses extracted by simple exponential fits to the diagonal 
components of the correlation function both for P = 1 and P = 3. Except for a possible 
small additive renormalisation there is no evidence for a qualitative difference between the 
two cases. There was no observable signal in the off diagonal components. It is possible 
that our lattices are too small (A^ = 1000 and A^ = 2000), but this we feel is unlikely since 
previous calculations for P = 1 with this model (ref. [9-10]) have indicated the onset of 
scaling behaviour on systems of this size. We thus conclude that the magnitude of any 
such coupling in the effective spin action is indeed very weak. 

It has been argued [6], that the coupling of c > 1 matter to two dimensional gravity 
creates instabilities in the worldsheet metric so to establish an order parameter for such a 
transition we should consider the effective gravitational partition function having traced 
over the spin variables. Associated with this is the question of what purely gravitational 
observables we should concentrate on. In our discrete model the geometry is determined 
by the local structure of the (p graph. This in turn is reflected in the size and frequency of 
its loops. A large loop corresponds to a highly coordinated vertex in the direct lattice and 
hence large scalar curvature. Thus, a non trivial quantity to focus on is the probability 
distribution Q (/) for the lengths of loops / in the graph. Moments of this distribution are 
then simply related to expectation values of moments of the scalar curvature density on 
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the direct lattice. Thus we are led to consider observables r (/3, P) defined by 

r^(AP)=J]Q(/)(l-//6)^ / = 3,...,iV (3) 

I 

One might expect that the large k moments are more sensitive to the occurrence of large 

(macroscopic) loops in the lattice. Fig. 1 shows a plot of the k = 2 moment describing the 

fiuctuations in the local scalar curvature density as a function of f3 for P = (pure gravity) , 

P — 1, P = 2 and P = 3 using A^ = 2000. On the figure we also mark the result for a, D = 1 

pure gaussian model on the same lattice. The strong coincidence of the gaussian model 

result with the two spin curve suggests a large universality at work, the local structure of 

the graph appears to be insensitive to all but the central charge of the matter and not the 

details of its fields or action. Clearly the addition of Ising spins enhances the fiuctuations 

in the local curvature peaking strongly at large P in the vicinity of the Ising critical point 

{Pc = 0.77332..). Furthermore the height of this peak scales linearly with P and hence the 

central charge. This is in qualitative agreement with the results of [11]. We have checked 

that this observation remains essentially the same if a different measure is adopted for the 

sum over graphs (e.g the conformal weight determined by adding a term — §X]j^'^(^i) ^o 

the action). However there is nothing in this plot signaling any pathology as P increases 

through two. These observations support a recent matrix model study by Brezin et al. 

[12]. 

One possible explanation for this stems from the observation that any c > 1 model is 
expected to have an infinite number of relevant operators. Any lattice representation of 
the theory would, in principle, have to include lattice versions of these operators, whose 
couplings would all have to be tuned to certain values for the system to become critical 
[13]. Since it is not clear how to construct these lattice operators (a generalised Jordan 
Wigner construction ?) and even less how to truncate the set to a manageable number 
whilst preserving the central charge, we chose to work with our simpler model. Since our 
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susceptibility data when P > 1 are still consistent with the scaling expected at a continuous 
phase transition, any continuum theory defined in its vicinity is still a candidate for a c > 1 
gravity coupled theory. 

For the case of pure gravity, the behaviour of Q (/) at large / is well known [3]. 

g(/)~e"""^ ao = 0.288.. (4) 

If we make the plausible assumption that this holds for P > 0, and consider the case of 
large k, it is easy to show that 

k\ f 1 ^ ^ 



Hence, if we divide the Ising moments by the corresponding pure gravity ones and consider 
the quantities 

Mk (/3, P) = In V^ (6) 



We expect 



M.(ftP)^Hn(^^l (7) 



Thus the parameters of the asymptotic distribution can be probed by fitting M to /c as 
both (3 and P are varied. 

In fig. 2, we plot M^^ {P = 0.7, P = 3) as a function of k for an A^ = 2000 node graph. 
Clearly we see good evidence for the expected linear behaviour, a least squares fit yielding 
a value of ln(ao/a3) = 0.18(1), x = 0.1. Using these fits we can compare the gradients 
In (ao/ctp) for P = 1, P = 2 and P = 3 from the A^ = 2000 node ensemble with P in the 
range 0.3 — 1.2 (figs. 3,4,5.). The P = 3 and P = 2 curves now show a prominent peak 
close to the critical coupling, indicating a rapid decrease in the associated ap, whilst the 
single Ising model exhibits a much weaker structure consistent with a broad plateau in 
the same region of coupling constant. Furthermore the peak height increases from 0.17(2) 
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to 0.18(2) for three Ising models on going from N = 1000 nodes to A^ = 2000, and from 
0.09(2) to 0.13(2) for P = 2. This is in contrast to the single spin case where the peak falls 
from 0.12(4) to 0.07(2) on going to the larger lattice. Whilst our errors are clearly large 
it is possible that our results are hinting at a qualitative difference in the thermodynamic 
behaviour of P < 1 and P > 1. Clearly high statistics runs on larger lattices will be needed 
to resolve this issue unambiguously. One possible scenario would be that the exponential 
behaviour in the tail of Q (/) for P < 1 gets replaced by a power law for P > 1. This would 
then ensure that sufficiently high moments of the scalar curvature would formally diverge 
signaling a breakdown of the model. However at present such ideas are merely speculation 
of how any possible continuum instability might manifest itself in the lattice models. 

It is perhaps surprising that the effect of adding Ising spins on the intrinsic geometry is 
rather similar to that of gaussian scalar fields. In the latter case the effective gravitational 
partition function is given by 

Z' = J]det"f (^-A^j (8) 



G 



N 



Here D counts the number of bosons (embedding dimension) and the scalar Laplacian 
A^- = Gij — qSij. Using this representation it has been shown that branched polymers (in 
the sense of the intrinsic geometry) dominate in the P> — ^ oo limit. 

In contrast, the critical region of the single Ising model is usually represented in terms 
of a Majorana fermion. On integration over this field one would have expected to find a 
positive power of an appropriate fermion determinant. Thus the dependence on graph G]\[ 
could be expected to be quite different from the bosonic case which is not what we see. 

We can gain some insight into this situation by considering the following (standard) 
representation of the Ising free energy F (/3, G) (eqn. 1) in terms of a continuous scalar 
field (pi- 

7 



If the dependence of F on graph G]\f could be evaluated exactly we could determine 
which graphs dominated the effective gravitational partition function in the Ising case. 
Unfortunately an exact solution is not possible and we may only compute F (/3, G) in some 
approximate way. The one loop approximation corresponds to expanding the field (p about 
some stationary configuration and integrating over small fluctuations. The stationary 
configuration is determined by the usual mean field equation 

G-^<j)j=l3tanh{<j),) (10) 

Although the generic graph G]\[ is random it still possesses a constant coordination number 
g = 3, so it is consistent to seek a homogeneous solution to this equation 

(j)^/l3q = tanh(i)^ (11) 

Now expanding the action about this stationary value (p to quadratic order and doing the 
resulting gaussian integrals leads to an expression for the free energy to one loop of the 
form (notice that the ill-defined det (G) piece has cancelled) 

F = Fi + F2 
Fi = N (in cosh (f)^ - ^ tanh^ cf)^ 

Notice that the mean field piece Fi does not depend on graph G. In the limit /9 -^ F2 
and Fi vanish and we recover the pure gravity result. Conversely, as /3 ^ 00, F2 gives 
no contribution whilst Fi contributes the {G independent) expected piece e^' '^i . For 
general (3 we can recast the F2 term as 

F2 = — Tr In [m^ - A^) (13) 



where A is the scalar Laplacian on the graph and m is given by 

2 cosh (j) 

Thus at /9 = /9c = I/q the mass m vanishes and a massless bosonic field emerges. 

In this situation the number of Ising models P plays a role analogous to the dimension 
D in the gaussian case and branched polymers dominate for large P. However, since we 
expect fiuctuation effects to be very important for two dimensions in determining the exact 
form of the effective action, the most we can hope for our gaussian approximation is that it 
indicates the sort of systematic trend we should expect. We feel that it is at least consistent 
with our numerical results. 

In conclusion, we have presented results from a simulation study of noninteracting 
Ising models on a dynamical random lattice. This was motivated by a desire to try to 
derive a better understanding of the coupling of c > 1 matter to two dimensional quantum 
gravity. Whilst we find no strong signal for a breakdown of the model in this region, 
we have perhaps pinpointed the relevant observables which might provide an unequivocal 
resolution of the question given a larger scale study. 

We have discussed some of the problems associated with a lattice formulation of the 
problem, and demonstrated within the context of a loop expansion an interpretation of at 
least some of the systematic trends we see in our Monte Carlo data. 
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FIGURE CAPTIONS 

[1] /c = 2 moment of the loop distribution for A^ = 2000 versus /3. P = 1 hy x, P = 2 
by o, and P = 3 by o. Solid lines indicate pure gravity (P = 0) result together with the 
D — 1 gaussian model result. 



[2] Mfc [13 = 0.7, P = 3) versus k, N = 2000. 
[3] In (^] versus (3, N = 2000. 
[4] In (^) versus P, N = 2000. 
[5] In (^) versus (3, N = 2000. 



p 


mp=i 


mp=3 


0.50 


0.520(4) 


0.530(1) 


0.60 


0.334(2) 


0.338(1) 


0.65 


0.245(2) 


0.247(1) 


0.70 


0.151(1) 


0.162(1) 


0.75 


0.069(1) 


0.089(1) 


0.80 


0.029(1) 


0.035(1) 



Table 1. Masses obtained from least square fits to simple exponentials, the fitting windows 
(which are the same for P = 1 and P = 3) being determined from the effective mass plots. 
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